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We discuss the separation of isospin-symmetric and isospin-breaking contributions in the hadronic 
observables within the framework of QCD plus QED. Further, we briefly review some recent work 
on the low-energy hadron phenomenology, in which the isospin-breaking effect plays a prominent 
role. 
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1. Introduction 

Strong interactions are described by the Lagrangian of pure QCD, whose only free parameters 
are the strong coupling constant and the quark masses. Most of theoretical predictions, which are 
made in the framework of pure QCD, assume, in addition, isospin symmetry, i.e. the masses of 
the up and down quarks are taken equal. In Nature, however, strong interactions do not come 
alone: the low-energy processes with the participation of hadrons contain electromagnetic effects 
- consequently, these processes should be described by the Lagrangian of QCD plus QED (for 
simplicity, we neglect weak interactions here). Therefore, in order to meaningfully compare the 
theory with the experiment, one should be able to "purify" experimental data with respect to the 
isospin-breaking corrections. This task is further complicated with the fact that the perturbation 
theory can not be applied to describe hadronic processes at low energies and one has to resort to 
the other techniques, be these the lattice QCD or Chiral Perturbation Theory (ChPT). 

Specifically, one faces the following problem. A physical process is described by the La- 
grangian of QCD plus QED, whose parameters are the strong coupling constant g, the electromag- 
netic coupling constant e and the (running) quark masses m u ,md,m s , ■ ■ ■. We define the isospin- 
symmetric world where e = and the masses of the up and down quarks are equal. The parameters 
of the isospin-symmetric world are the strong coupling constant g and the quark masses m u ,m s ,---, 
with m u = thd (the bar over a symbol denotes the isospin-symmetric world). In general, the follow- 
ing questions should be answered: 

i) How are the parameters g,e and m u ,md,m s , ■ ■ -, on the one side, and g and m u ,m s ,-- -, on the 
other side, related? 

ii) How does the separation of isospin-symmetric and isospin-breaking effects translate to the 
level of the effective chiral Lagrangian, which is used to describe hadronic processes at low 
energy? 

iii) Once the prescription for separation of the isospin-symmetric and isospin-breaking effects is 
set, what is the systematic procedure of purifying the hadronic observables with respect to 
the isospin-breaking corrections? 

Recently, some aspects of the problem have been addressed, e.g., in Refs. [1-4]. Below, we 
mainly follow the line, adopted in Ref. [3]. Instead of a general reasoning, for illustrative reasons, 
we concentrate on a particular example. 

2. "Purely strong" pion decay constant 

In this section, we shall study the question, how to extract the "purely strong" pion decay 
constant F n from experimental data on the decay of a charged pion with or 1 real photon in the 
final state. This problem has a long history, starting in the pre-ChPT era (see, e.g. [5]). In ChPT, 
the decay width of this process has been calculated at one [6] and two [7] loops. Below, we give an 
expression of the decay rate at one loop, evaluated in ChPT with three flavors [6] 
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where ni( denotes the lepton mass, M n + ,M n a,M K + ,M K o are physical pion and kaon masses, and 
M 2 = 2Bom, M\ = Bo(m s + m) with m = A (m l( + m,/). Further, Fo denotes the pion decay constant 
in the chiral limit, Bo is proportional to the quark condensate in the chiral limit, L\ denote the strong 
low-energy constants (LECs), E r is a certain linear combination of the electromagnetic LECs, the 
function H{z) stands for the contribution of the photon loops, and pt is the scale of the dimensional 
regularization. Finally, Gp and V u d stand for the Fermi constant and the element of the Cobayashi- 
Maskawa matrix. 

The reasoning goes as follows. In the isospin-symmetric world, the quantity Fq can be related 
to the charged pion decay constant F n through the well-known expression [8] 

F% = F °{ 1 + f\ {L * {Ml + 2M * } + LWl) " 32^ ln F + ^ W) } ' (22) 

Here one has implicitly assumed that passing to the isospin limit amounts to setting e = 0, = 
m u = m, whereas all other parameters of the theory stay put. Further, substituting this expression 



into Eq. (|2. 1|), it is seen that one may extract the exact value of F n from the measured value of the 
decay rate, provided one makes reliable estimates for the electromagnetic LECs contained in E r . 

In order to demonstrate that this procedure is ambiguous in general, we calculate the pion 
decay constant in a model where, unlike QCD, explicit calculations can be performed. To this end, 
we invoke the linear a -model to one loop. In the limit when the mass of a becomes much larger 
than the pion mass, the low-energy structure of the theory can be described by ChPT - with the 
specific values of the LECs which are determined by the underlying Lagrangian of the a-model. 

In the following, we mainly follow Ref. [3]. The Lagrangian of the linear a-model with 
electromagnetic interactions is given by «5f = «5fo + -^ct, where 

J?o = \d^ T d^ + ^f^- 8 -^ T ^) 2 + c^ + ^f{Q^ T (Q^) 

- ^mfiQm 7 *)- 1 ^^-^^) 2 , (2.3) 

where = denotes the 4-component spin-0 field, is the electromagnetic field, F^ v is 

the electromagnetic field tensor, ^ is the gauge parameter and Q stands for the 4x4 charge matrix, 
whose non-zero components are Q \i = —Q21 = 1- Further, the co variant derivative of the field (p 
is defined as d^ip = d^ty + {F 11 +eQA fl )(f), where external vector and axial- vector fields are given 
by = a'^ , F'^ J = —e'^v 1 ^. Finally, the counterterm Lagrangian Jz? ct = £^ =1 j8 fl ^ a includes all 
operator structures, which are necessary to remove divergences to one loop (for more details, see 
Ref. [3]). The isospin-breaking parameters parameters 8m 2 , 8g are counted at order e 2 . 

The spontaneous chiral symmetry breaking in the model occurs, when m 2 > 0. In this case, 
the parameter c, which describes the explicit breaking of chiral symmetry, determines the mass of 
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F n=Z n{ WW + ^\aaa, } 



Figure 1: The neutral pion decay constant at one loop in the linear c-model. 

the neutral pion (the charged pion mass receives an additional piece proportional to e 2 , see [3]). 
The mass of the a-meson, M a , does not vanish in the chiral limit. 

In order to evaluate the pion decay constant, we calculate the matrix element of an axial current 
between the vacuum and the neutral pion state to one loop. This matrix element is determined by 
the diagrams shown in Fig. [l| At the next step, we perform the limit M a — > in the resulting 
expressions, and compare the answer with the neutral pion decay constant, calculated in the two- 
flavor ChPT [9] 

F n o-F< 1- In — y + -=~~ 14 — e > ak t } , (2.4) 



l6n 2 F 2 pL 2 F 2 



where F, l\ are the ,S£/(2)-analogs of Fq,L\, and c\k\ stands for a some linear combination of the 
electromagnetic LECs. The matching gives 

Now, we must move one step further and relate the quantity F to its "purely strong" counterpart, 
defined in the limit e = 0. This procedure is straightforward in Quantum Mechanics but becomes 
obscure in the framework of the Quantum Field Theory, where one has to deal with the renormal- 
ized parameters. The renormalization group (RG) equations, which govern the running of these 
parameters, to one loop are given by (see Ref. [3]) 

dm 2 1 

^luT = W- <(3s + Ss) +sSm ~> {2 - 6) 



Whereas the RG equations in the isospin-symmetric world are obtained by setting e = in 



Eq. (2.6) 



dm 2 




IX 

d\l 








djJL 


= — f 



(2.8) 
(2.9) 

The crucial point here is that the procedure of switching off the electromagnetic corrections unam- 
biguously predicts the running of the parameters with respect to the renormalization scale }i, but 
not the value of these parameters at a given fi. The parameters g(ju), m 2 {\i) can be, in principle, 
fixed from the experimental data. However, this is impossible for g(jJ.), m 2 (jj.), since there are no 
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Figure 2: Matching of the strong coupling constant in full theory and in case when the electromagnetic 
interactions are switched off. The figure illustrates the uncertainty related to the choice of the matching 
scale. 



data with electromagnetic interactions switched off. For this, in order to fix the values of g(p.), 
m 2 (/i), one has to invoke additional conventions, rendering the definition of the isospin-symmetric 
limit convention-dependent. For example, one may impose boundary condition on the solutions of 



Eq. (2.8) at some jj. = which will be hereafter referred to as the matching scale. Most easily, 



this can be done by requiring g(}X\ ) = g(}X\ ) and m 2 {}Xi )=m 2 (ni). Then, from Eqs. ( p76[ ) and ( p78| ) 
one obtains 

m 2 Gu) = m 2 Gu;Ati) + iln — (8gm 2 + g8m 2 ), (2.10) 
47T Z Hi 

8 (H) =g(n;fii) + ^\n^- g 8g. (2.11) 

Here, we explicitly indicate the /ii -dependence of the parameters of the theory in the isospin sym- 
metry limit. This dependence is illustrated in Fig. || where it is shown that, changing the matching 
scale from }X\ to jj, 2 leads to the change of the value of g(jj,) at a fixed }X. 



From Eq. (2.5) one now obtains (remember that the quantity F , by definition, refers to the 



isospin-symmetric world) 



m 

7I\ 



li-js^+jsy (212) 

\ 167T 2 jU 2 327T 2 J 



One may directly check that = 0, as it should. However, the quantity F depends on the 
choice of the matching scale since ;Ui J~ ^ 0. This is the ambiguity inherent to the definition 
of the isospin-symmeUy limit. There is no way to avoid this ambiguity. On the other hand, a 



particular combination of F and the electromagnetic LECs, which is present in Eq. ( |2.5D , is free of 
this ambiguity 1 . 



'This way of reasoning applies also to the quark masses which are present in the Lagrangian of ChPT: these are 
the parameters of pure QCD and not of QCD plus QED. For instance, u— and d— quark masses in the Lagrangian of 
ChPT have the same RG running, so that their ratio is RG invariant. In analogy to F , these masses are also convention- 
dependent. 
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It turns out that within the framework of the linear a-model it is possible to make a numerical 
estimate of the size of this uncertainty. This happens because the quantities 8g and 8m 2 , which 



determine the /ii dependence of the parameters in Eq. ( |2.10[ ), also enter the expression of the 
charged and neutral pion mass difference. Using this fact, one obtains the following estimate 
within the linear a-model [3] 

F(ni = 1 GeV)-F(/ii =0.5 GeV) -0.1 MeV. (2.13) 

That shows that the effect is not purely academic (for example, one may compare this with the 
uncertainty in F n , quoted recently in Ref. [10]). 
The following conclusions can be drawn: 

i) The separation of the electromagnetic and strong interactions is not an unambiguous oper- 
ation in Quantum Field Theory. We have seen this explicitly on the example of the linear 
a-model at one loop, and we have no reason to believe that this will be different in QCD. 

ii) As explained above, this ambiguity can be fixed, e.g., by setting boundary conditions in the 
RG equation. As a result, the parameters of the isospin-symmetric theory will depend on the 
choice if the matching scale /ii . Alternatively, one may use observable quantities for fixing 
the ambiguity. In the example considered in this section, it suffices to set the pion decay 
constant and the pion and the a masses in the isospin-symmetric world to a given values 
(known as reference values). This determines the quantities g, m 2 and c uniquely. All other 
observables are expressed through these parameters. 

iii) Returning to the issue of extracting the value of F % from the data, we see that the extraction 
with an arbitrarily high precision is not possible, because the quantity F (and hence F n which 
coincides with F in the chiral limit) depends on the matching scale jJ.\ . 

iv) Using ChPT with virtual photons, which was originally introduced in Refs. [11, 12], to cal- 
culate isospin-breaking corrections to the "purely hadronic" observables implies that the 
hadronic LECs stay put in the isospin symmetry limit. In the linear a-model this means that 
the LECs are expressed through g and m 2 , so that they depend on ;Ui. On the other hand, the 
physical quantities in the real world do not depend on jX\. Consequently, the jii -dependence 
in the strong LECs must be compensated by a redefinition of the electromagnetics LECs. In 
the example with extracting F % from the data this would mean that there exists an inherent 
uncertainty which does not allow us to fix the electromagnetic LECs - and hence F n - with 
an arbitrary precision. 



3. Isospin-breaking effects in hadronic observables 

From the previous section one may conclude that the definition of the (hypothetic) isospin- 
symmetric world, which sets the reference point for the calculation of the isospin-breaking cor- 
rections, is ambiguous. On the other hand, there exist observable isospin-breaking effects, where 
the above ambiguity does not matter. Below, we want to briefly review some recent work where 
the isospin-breaking effects play a prominent role, and highlight the issues which are related to the 
discussion given in the previous section. 
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3.1 Cusps in the K -^3n decays 

Recently, NA48/2 collaboration at CERN observed a pronounced Wigner cusp [13] in the in- 
variant mass distribution for the three-pion decays of charged kaons [14]. Cabibbo in Ref. [15] 
proposed that measuring this cusp, which is due to the presence of the charged and neutral pion 
mass difference, can be used for extracting the values of the nn scattering lengths (originally, the 
cusp has been predicted in Ref. [16], see also [17] where the cusp in the %% scattering amplitude is 
discussed). However, the accuracy of the one-loop representation of the K^3% decay amplitudes, 
which is given in Ref. [15], does not suffice to fit the high-precision experimental data. To this end, 
in Ref. [18] one has used analyticity and unitarity of the 5-matrix to obtain a theoretical representa- 
tion of the charged and neutral kaon decay amplitudes, valid up to and including second order in the 
%% scattering lengths. The authors of Ref. [19] have merged the above approach with ChPT, which 
is only used to calculate the real part of the K — > 3n decay amplitude at one loop. In Ref. [20], the 
charged pion decays have been studied to the same accuracy within the non-relativistic effective 
Lagrangian approach, which turns out to be a most systematic and convenient tool to address the 
problem. In particular, it automatically includes the strictures imposed by analyticity and unitarity. 
Furthermore, the inclusion of the electromagnetic effects proceeds straightforwardly within this 
approach [21]. Using two-loop representation of the decay amplitudes, a precise determination of 
the % % scattering lengths is possible [22]. Finally, we mention two more approaches to the same 
problem. Ref. [23] uses the dispersive approach. Recently, there have been attempts [24] to ad- 
dress the problem in a framework which may be considered as a merger between ^-matrix theory 
and a conventional quantum-mechanical framework used to study Coulomb interactions in the final 
state. In general, this framework is contradictory and incomplete, since, e.g., it fails to reproduce 
the correct two-loop structure of the decay amplitude, as well as the full set of isospin-breaking 
corrections to it. 

It should be pointed out that the non-relativistic approach can be straightforwardly general- 
ized to study other three-particle decays, like the decays of the neutral kaons [25], (for the ex- 
periments on neutral kaons, see Refs. [14,26]), as well as tj — > 3n decays [27] and r\' — ► r\2n 
decays [28]. From the recent developments, we also highlight a comprehensive calculation of the 
isospin-breaking corrections at order e 2 (m^ — m u ) in r\ — ► 371 decays, carried out in ChPT to one 
loop [29]. These calculations complement earlier calculations at order e 2 m, carried out in Ref. [30]. 
For more theoretical and experimental developments on the cusps, see, e.g., contributions to this 
workshop [22, 31]. 

We would like to stress that the representation of the decay amplitude, which is considered 
in the present subsection, concerns the real world and the observable isospin-breaking effects (like 
the emerging cusp) only. No reference is made to the idealized world with no isospin breaking 
until the very end, when the relation between the physical %% scattering amplitudes in different 
channels and the 5-wave %% scattering lengths ao,a2 is established (see Ref. [20] for more detail). 
Fortunately, the isospin-breaking corrections, contained in this relation, at leading order can be 
expressed by the charged and neutral pion mass difference and are thus parameter-free [20]. The 
electromagnetic LECs, which are the source of ambiguity discussed in the previous section, appear 
only at next-to-leading order. Their contribution are expected to be small and well-controlled. 

Finally, we would like to emphasize that ChPT [32], combined with Roy equations, allows 
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one to make very precise predictions for the values of the S-wave scattering lengths in elastic %% 
scattering [33]. For this reason, confronting the experimental results with these predictions, one 
may extract important information about the fundamental properties of QCD at low energy. For 
instance, the above predictions for the scattering lengths have been made in the assumption that 
the spontaneous chiral symmetry breaking in QCD proceeds according to the standard scenario 
(with the large quark condensate). If the experimentally measured scattering lengths significantly 
deviate from the theoretically predicted values, this would indicate that the symmetry breaking in 
QCD follows a different path. 

3.2 Isospin-breaking corrections is K e 4 decays 

K e 4 decays represent an important source for the determination of the %% scattering phase 
shift and scattering lengths [22, 34, 35]. The crucial observation which allows one to link the K e 4 
decay amplitude to the %% elastic scattering amplitude is that, according to the Watson theorem, 
the phases of both amplitudes coincide. Watson theorem, however, assumes isospin symmetry, 
which is violated in Nature. It is therefore justified to ask, how large are the corrections in the %% 
scattering phase, which is extracted from the K e 4 decays. 

In fact, in the analysis of the experimental data one has already removed part of these correc- 
tions, applying the Coulomb factor and using the program PHOTOS [36] - see Ref. [35] for details. 
This procedure definitely leaves out some of the isospin-breaking corrections, e.g., those caused by 
the mass differences in the isospin multiplets, or direct effects related to the quark mass difference 
Wd — m u- Thus, one is faced with an alternative: either one chooses to deal with the "raw" exper- 
imental data and removes a full set of isospin-breaking corrections which can be systematically 
calculated in ChPT, or one uses the data already "corrected" by PHOTOS and removes additional 
isospin-breaking effects, assuming that this can be done separately. Although the first option is def- 
initely more appealing, it requires a major calculational effort (see Ref. [37] for an earlier work in 
this direction). Recently, the calculation of the isospin-breaking corrections, using the "corrected" 
data, has been carried out in Ref. [38] (see also Refs. [39] which attempt to calculate the effects 
caused by the pion mass splitting and electromagnetic corrections). The findings of Ref. [38] can 
be summarized as follows: 

i) The isospin-breaking correction have been calculated to one loop in ChPT. Virtual photons 
have been neglected. The charged pion mass is used as the reference mass for the isospin 
symmetric world. Both the corrections coming from the pion mass difference and the effect 
proportional to — m u , which originates from the n° — f] mixing phenomenon, are taken 
into account. Both these effects are of a comparable size. 

ii) Watson theorem holds only in case of isospin symmetry. When isospin symmetry is violated, 
the phases in the K e 4 decay amplitude and in the nn elastic amplitude differ. Moreover, they 
are not related to each other: for example, the decay amplitude contains some electromag- 
netic LECs which are not present in the %% amplitude. 

iii) If isospin symmetry is violated, the phase of the K e 4 decay amplitude has a cusp below 
threshold. The phase no more vanishes at threshold s = 4M^ + . 
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Figure 3: The difference between the phase in the K e 4 amplitude y/b and the nn elastic S-wave scattering 
phase with full isospin / = in the isospin limit 5q. The blue band shows the error caused by the uncertainty 
in the parameters of the calculation. M nK denotes the invariant mass of the charged pion pair. 



iv) It has been shown that the isospin breaking effect is quite large (see Fig. ||). Moreover, these 
corrections strongly affect the values of the S-wave %% scattering lengths, which are ex- 
tracted from the experimental data [22]. Only after the isospin breaking is taken into account, 
the experimental output for the scattering lengths agrees with the theoretical prediction [33]. 



3.3 Isospin breaking in the kN scattering and in the pionic hydrogen observables 

The experiments on hadronic atoms [4(M-2] provide a beautiful possibility to directly extract 
the hadron-hadron scattering lengths from the experimental data without using the extrapolation to 
the threshold. However, the accuracy of the extraction of the scattering lengths critically depends 
on the possibility to gain a control on the isospin breaking corrections in the atomic observables. 

In Ref. [43] the relation of the ground-state energy level shift in the pionic hydrogen and the 
elastic %~ p threshold scattering amplitude is obtained at next-to-leading order in isospin break- 
ing. Further, in Ref. [44] a similar relation between the width of the ground state and the charge- 
exchange 7T /? — > n°n threshold amplitude has been derived. The accuracy of these relations com- 
fortably matches the existing experimental precision. However, the question of calculating the 
isospin-breaking corrections for the threshold amplitudes, in order to relate these quantities to the 
S-wave TtN scattering lengths, has proven to be more difficult. One invokes ChPT to this end. 
In Ref. [43], the isospin-breaking corrections to the %~ p amplitude has been obtained at 0(p 2 ), 
using the charged pion and proton masses as reference masses in the isospin-symmetric world. In 
Ref. [45] these calculations have been extended to 0{p i ). Finally, in Ref. [44] the isospin-breaking 
corrections to the charge-exchange amplitude have been evaluated at o(p 2 ). For the recent review 
on the subject, see, e.g., [46,47]. 

There are, generally, two problems, related to the calculation of these corrections: 

i) Using leading order values for the corrections does not always suffice. Loop corrections can 
be sizable, see, e.g., Ref. [45]. 
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ii) The lowest-order LECs (most notably, the electromagnetic LECs) are poorly known. En- 
tering the result already at leading order, these LECs are responsible for the bulk of the 
uncertainty in the calculations. 

Although the calculations at higher orders in ChPT can not improve on item ii), they can 
definitely check the convergence of the chiral expansion and lead to a more reliable estimate of 
the isospin breaking corrections, see item i). This is the major justification for the comprehensive 
study of the isospin-breaking effect in all physical nN scattering amplitudes at threshold, which 
has been undertaken in Ref. [48] at 0{p i ) in the relativistic baryon ChPT. A careful analysis given 
in Ref. [48], which involved both extensive analytic calculations and a significant numerical effort, 
will allow one to carry out the extraction of the nN scattering lengths from the pionic hydrogen 
data in a more controllable manner. 
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